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ABSTRACT 

We calculate the proton distribution function due to heating by subsonic (Alfvenic) 
turbulence in a weakly magnetized collisionless plasma. The distribution function is 
nonthermal. For non-relativistic energies, it is an exponential of the magnitude of the 
proton velocity. For ultra-relativistic energies, it can be characterized as a power law 
with a momentum-dependent slope. 


Subject headings: plasmas — MHD — turbulence 


1. Introduction 

The proton distribution function (PDF) in a collisionless plasma depends on the history of the 
plasma. If, for example, the plasma is heated by shocks, the PDF is believed to be a power law at 
high energies, while for the thermal particles, the form of the PDF is not thought to be universal 
(e.g., Blandford & Eichler 1987). In this paper we calculate the PDF for a weakly magnetized 
plasma (the Alfven speed va is less than the ion thermal speed Vi) heated by subsonic (Alfvenic) 
turbulence. 

Although we do not have a fully quantitative theory for the turbulence (§2), we are still able 
to calculate the distribution function (§3). This is due to a special feature of particle heating by 
low frequency, subsonic, waves (§3.1). In §4 we mention a possible application of our results to 
accreting black holes. We summarize our results in §5. 


2. Proton Heating by Subsonic Turbulence 

The nature of subsonic turbulence in a collisionless weakly magnetized plasma is not fully 
understood. A plausible scenario, based on Goldreich and Sridhar (1995), is as follows. Large scale 
motions (larger than the ion Larmor radius p) are nearly dissipationless and incompressible. If 
the characteristic velocity of the energy containing eddies (eddies of the largest scale L) is small, 
V <; va <C Vi, the dynamics is adequately described by reduced MHD (incompressible MHD without 
slow modes, e.g. Strauss 1976, Ng & Bhattacharjee 1996). As in hydrodynamics, the turbulent 
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energy cascades to small scales where it is dissipated. The cascade is described by a Kolmogorov 
5/3 law, but it is strongly anisotropic, with most of the energy residing in quasi-perpendicular 
perturbations. Typical wavevectors k are nearly perpendicular to the local magnetic held, so that 
perturbations are strongly elongated along the local magnetic held. The energy per unit mass 
in eddies of perpendicular wavenumber ~ k± is e 2 / 3 fc_|_ 2//3 , where 6 r '- / V 3 /L is the net energy 
dissipation rate per unit mass. The degree of anisotropy is determined by the requirement that the 
Alfven frequency not exceed the nonlinear frequency, which yields kn <; fe^ 3 L -1 / 3 . 

“Viscous” heating occurs when the perpendicular size of the eddy is ~ p (Quataert 1998, 
Gruzinov 1998, Quataert & Gruzinov 1998). The small-scale eddies are essentially Alfven waves. 
Alfven waves with a perpendicular wavenumber ~ p~ l have a non-zero parallel magnetic held 
perturbation (in contrast to long wavelength perturbations). As a result, particles interact with 
waves through “magnetic mirror” forces, i.e., through the coupling of the magnetic moment of the 
particle and the parallel gradient of the magnetic held. A charged particle moving along the held 
will be randomly accelerated and decelerated along the held, and will therefore diffuse in parallel 
velocity space. 

Strong parallel diffusion requires |vy | <; va, where vn is the component of the particle’s velocity 
along the magnetic held. This can be understood as follows. A parallel magnetic held perturbation 
of frequency c o and parallel wavenumber kn will accelerate particles that satisfy the resonance 
condition uii = c o/k\\. If the turbulence were strictly a superposition of linear Alfven waves, the 
resonance condition would imply that only particles with |u||| = va are accelerated (because Alfven 
waves have the dispersion relation oj = |A:ii|u^). 

It is likely, however, that the strong turbulence of interest is not fully describable as linear 
Alfven waves. In this case, at the smallest scales ( k± ~ p _1 ), perturbations of all frequencies 
u £ u n i and wavenumbers kn ^ c v n i/vA will exist, where oj n i ~ e 1 / 3 ^ -2 / 3 is the nonlinear frequency 
on scales ~ p. The quantity u:/k\\ then takes on values between 0 and oo. This does not, however, 
mean that a particle with an arbitrary parallel velocity will be accelerated. Acceleration of high 
parallel-velocity particles, |uy| > va , will be insignificant because perturbations resonating with a 
high parallel-velocity particle must have small kn. These perturbations will exert negligibly small 
magnetic mirror forces for two reasons. The parallel force is proportional to the parallel gradient 
of the parallel magnetic field perturbation (see eq. Q). In addition, the parallel magnetic field 
perturbation is itself proportional to the parallel gradient of the electric field perturbation, because 
an incompressible two dimensional plasma flow with fey = 0 does not perturb the parallel magnetic 
field.[] Magnetic mirror forces are therefore ~ k?r, perturbations with small fey are consequently 
negligible. 


1 Formally, this follows from the small fey limit of the plasma dielectric tensor. 
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3. The Evolution of the Proton Distribution Function 


Diffusion in parallel velocity due to turbulent heating can be described by the equation 
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where Dy is the parallel momentum diffusivity. Anisotropies in momentum space are quickly 
erased by high-frequency small-scale electromagnetic instabilities. Parallel momentum diffusion is 
then equivalent to an isotropic diffusion in momentum space, 
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where 

D = \ f i d W 2 D\\ (3) 

and p is the cosine of the particle’s pitch angle (the angle between the velocity vector and the local 
magnetic held). 


Suppose that protons are strongly heated, so that their final thermal energy is much larger 
than their initial thermal energy. In this case, the final distribution function depends primarily 
on the properties of the diffusion under consideration, but only weakly on the initial distribution 
function. If D does not depend on p , the final distribution function will be Maxwellian.^] As we 
show below, for our problem D oc qp, where 7 is the particle’s Lorentz factor. Since D is not 
constant, diffusion establishes a non-Maxwellian PDF. 


3.1. The Diffusion Coefficient 


We first calculate the parallel diffusivity, D y. The parallel motion of a particle interacting with 
low frequency perturbations is described by 


dv \| 
dt 


Vy B 

~2 IT' 


( 4 ) 


where we have assumed 


<C v±. Since p» = 'ymv n, the parallel momentum diffusivity is 

D\\ oc ^ 2 v\^{v\\/va)- (5) 


<P(x) is a positive dimensionless function, which is ~ 1 for |x| <; 1 , and rapidly approaches zero at 
larger |x|. This is because only particles with sub-Alfvenic parallel velocities efficiently interact with 
the turbulence (see § 2 ). d>(x) cannot be calculated in detail because the spectrum of turbulence at 


2 The Green’s function of equation (^) with D = const, is a Gaussian. 
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small scales, where nonlinear and kinetic effects are both important, is unknown^] Because we are 
interested in va Vi, however, our hnal answer will be independent of <J>. 

The isotropic diffusivity, D, can be obtained from D\\ using equation (|), which yields 

Dcz'y 2 v 4 J dpp 2 (l — p 2 ) 2 ^{pv/ va) oc 7 2 f. (6) 


3.2. The Distribution Function 


For non-relativistic protons, D oc p. Solving equation (|2j) for an initial PDF fo(p) oc 5(p) gives 

f(p) = Jrp— c 3 /9 ex P(~2p/ VTrn) , (7) 

7r(J m) 6 i z 

where the “temperature” is defined by the usual expression < p 2, >= 3 Tm. Although the PDF 
is nonthermal, the number of suprathermal protons is exponentially small. Subsonic turbulence is 
therefore intrinsically inefficient at accelerating particles in the non-relativistic limit. 


At ultra-relativistic energies, D ocp 2 . Solving equation 

-a(p) 


for f 0 (p) oc 5(p - pi), gives 


where 


f(p) oc — 
\Pf 


a{p) = I + lQ g (P/P/) 

U 2 + log ( Pf / Pi ) 


( 8 ) 


(9) 


and pf is the mean momentum of the final PDF. At relativistic energies, there is a significant 
population of particles above the mean energy of the plasma. Realizing this efficient acceleration, 
however, requires that the bulk of the plasma be heated to at least mildly relativistic energies, 
which limits its practical significance. 


Using equation (]2]) and equation (|6|), we numerically calculated the distribution function for 
a plasma that was heated from non-relativistic to mildly relativistic energies. Figure 1 shows the 
calculated distribution in energy space, dN/dE, where E = (7 — l)mc 2 . Also shown is a Maxwellian 
with a temperature of 0.2GeV, which has the same mean energy per particle. The viscously heated 
PDF has a noticeable nonthermal tail of relativistic particles. 


4. Application to Accretion Flows 

A possible application of our results is to hot, quasi-spherical, collisionless plasmas which have 
been proposed to exist near accreting black holes. By spherical accretion we mean any of the 

3 For linear MHD Alfven waves, $(*) oc <5(|ar| — 1). 
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proposed accretion models, Bondi (Bondi 1952), ion tori (Rees et al 1982) or advection-dominated 
accretion flows (Narayan & Yi 1995). Our calculated PDF may be relevant to these gravitationally 
confined plasmas (which must have va ^ D and, as we have argued elsewhere, are likely to have 
va Vi). For protons with roughly virial energies near the Schwarzschild radius of a black hole, the 
PDF calculated in the previous section is noticeably non-Maxwellian, with an excess of high-energy 
protons. 

In such hot plasmas, even a Maxwellian PDF leads to interesting levels of gamma-ray emission 
due to pion production in proton-proton collisions (Mahadevan et al 1998; Ginzburg & Syrovatskii 
1964 give a clear description of how to estimate the rates). Our non-Maxwellian PDF should 
produce distinct observable features in the gamma-ray spectrum. This will be calculated in a 
separate paper since theoretically estimating the PDF for the purposes of gamma-ray emission 
entails (in addition to the present analysis) a number of additional complications, e.g., adiabatic 
heating of particles as well as the possibility that the flow may capture and compress ambient 
cosmic rays (Meszaros 1975) or pre-accelerated protons from a Bondi-Hoyle bow shock. 


5. Summary 

We have calculated the proton distribution function due to heating by subsonic (Alfvenic) 
turbulence in a weakly magnetized collisionless plasma. For non-relativistic energies, the PDF 
is nonthermal only to the extent that it cuts off exponentially, rather than as a Gaussian, with 
particle momenta. By contrast, at mildly relativistic or relativistic energies, a potentially observable 
“suprathermal tail” develops. 

Our results seem to be relatively robust. For plasmas heated to well above their initial thermal 
energy, the distribution function is determined solely by the momentum space diffusivity. The 
diffusivity we have calculated (eq. [j(|) should be valid on quite general grounds, provided the 
protons are indeed predominantly heated by low-frequency, subsonic, waves. 
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